Spiral structures in helical magnets 
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Structure and properties of two-dimensional spiral textures in helical ferromagnets have been 
studied. In these magnetic mediums have been predicted new types of periodical structures - spiral 
lattices. 



Different spin structures in magnetic systems are 
the matter of intensive theoretical and experimen- 
tal researches. Last years a study of spin textures 
in non-centrosymmetric helicoidal ferromagnetic crys- 
tals have aroused doubtless interest. Here, owing 
to Dzyaloshinskii-Moriya interaction, the ground state 
is one-dimensionally modulated phase - helicoid [1]. 
Thereby such spin ordering with edge magnetic dislo- 
cations was discovered in cubic crystal Feo.5Coo.5Si by 
methods of transmission electron microscopy in the pa- 
per [2]. In the following paper [3] the same methods 
were used for observing spiral (Swiss-roll-like vortex) in 
helimagnet FeGe with continual heating of the sample 
from 85K to 200K. Finally, skyrmion lattice, which was 
predicted many years ago [5, 6], was found recently [4]. 
To the present date there is no theoretical and numeri- 
cal analysis of two-dimensional structures in helimagnets 
except skyrmions and their lattices. 

In this paper by using analytical and numerical meth- 
ods for two-dimensional model we research the structure 
and properties of spiral textures, which were recently 
discovered experimentally [3]. An analytic formula for 
the structure of spirals outside the corn was calculated. 
It matches with numerical calculations and it gave an 
opportunity to research their properties including be- 
haviour in the magnetic field. Existence of single struc- 
tures in condensed matters (dislocations, disclinations, 
vortexes, skyrmions, etc) is as the rule accompanied by 
the formation of their periodical structures. We are pre- 
dicting new types of periodical structures in helical fer- 
romagnets - the Archimedean spiral lattices. 

The energy of the helical ferromagnetic cubic crys- 
tals (such as FeGe, MnSi [7]) includes contributions from 
the exchange interaction, external magnetic field H z and 
Dzyaloshinskii-Moriya interaction: 
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where n = (n x , n yi n z ) is a unit vector of magnetization: 
n = (sin6cos<j), sin6sin<j),cos6), n 2 = 1, (2) 




Figure 1: Calculated Archimedean spirals surrounded 
by a different structures. (a),(b) AS surrounded by a 
labyrinth structure, (c) AS surrounded by a skyrmion 
lattice in zero field; (d) - in the filed h = 0.5. 



M - saturation magnetization, a - exchange constant, D 
- Dzyaloshinskii constant. 

The Euler-Lagrange equations for energy functional 
(1) in polar coordinate system (r,ip): 
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where A is the Laplace operator and parameters 

h = H z /H D , H D =D 2 M/a, L D = 2ira/D (5) 

are the same as in the paper [8]. 

The method of energy functional (1) minimization was 
used for solving equations (3)- (4). Calculations were 
made on the high density grid: 2048x1024. So that 
the angle between the two directions of a vector in the 
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adjacent nodes did not exceed 10°. The nonlinear con- 
jugate gradient method with additive penalty functions, 
which are required for vector norm fixation (2), was used. 
To control the result accuracy the first and the second 
derivatives were calculated and then directly discrepancy 
in the equations (3)- (4). Algorithm was realized on the 
Nvidia CUDA architecture using video cards micropro- 
cessors (GPU) for massive parallel calculations and vi- 
sualization in real time. The usage of such technolo- 
gies allow to increase sufficiently the speed/accuracy of a 
number of numerical problems, including micro-magnetic 
ones [9]. As a result different types of spiral structures 
were found. They will be described further. Moreover, 
the analysis of obtained numerical results gave us an op- 
portunity to make analytical research of their structure 
and properties too. 

We found asymptotic spiral structure of solutions of 
(3)-(4). Thus for h = 0: 
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where N £ Z. The origin of logarithmic term is tied 
to the exchange spirals [10]. In the case D ^ it is 
necessary for minimization (1) to set arbitrary constant 
C2 = 0. Then solutions for n z = cos (9) are 2N-spiral 
domains, which are divided by Archimedean spirals. It 
is possible to show, that when h > asymptotic solution 
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where am - Jacobi amplitude, F - elliptic integral of the 
first kind, K = K(k) - complete elliptic integral of the 
first kind. Elliptic modulus h depends only on h and 
determines by equation 
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where E = E(x) complete elliptic integral of the second 
kind. It is noteworthy that the period of spiral turns 



L 



D 



:K 



4EK 



(9) 



and coincides with the period of helicoid [1, 11, 12]. 

If TV = 0, then solution conforms to the extreme case 
of kir -vortices [13]. For \N\ = 1, level lines n z = const 
are the Archimedean spiral (AS) (if \N\ > 1, then spiral 
is multiarm). 

By integrating the functional (1) for the leading terms 
(7), we find the energy, which is contained in k turns 
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where Eo = aM 2 , e c - dimensionless quantity of the core 
energy, which depends only on h and N. Core structure 
can not be determined by asymptotic method and energy 
e c was obtained by numerical method. 

It follows from (12) that the direction of AS curling 
and a number of arms does not influence spiral energy 
density at a distance from the core. As a result of this 
we studied only the simplest case N = 1. Such kind of 
AS was observed in the experiment [3]. 

As we can see from the numerical calculation, spi- 
ral textures as metastable states can exist surrounded 
by labyrinth structures (fig.la,b) or surrounded by 
skyrmions (fig.lc,d), like static spiral domains (SD), sur- 
rounded by labyrinth domain structure [14] or by bubble 
lattice [15, 16] in garnet ferrite films. In these films SDs 
may contain up to 30 turns [17] and being observed only 
as dynamical structures in some samples [18]. Magnetic 
dipole interactions play the key role in forming static SDs 
[19], but for the Archimedean spiral in chiral magnets the 
same as for other solitons with typical sizes about Lp, 
taking into account the influence of these long-ranging 
forces is not so critical [20]. 

On the fig. 1 we can see how does the AS changes its 
form a little under the influence of the surrounding. Sim- 
ilar spirals for the director field can be observed in the 
liquid crystals films [21-23]. However, director field is 
highly inhomogeneous over the film thickness and con- 
tains topological defects [24]. AS vector field is contin- 
uous everywhere and does not contain any singularities 
but it is topologically nontrivial. Localized AS (fig.lc,d) 
can be smoothly deformed into Belavin-Polyakov vortex 
[25] with unit topological charge. For instance, if the field 
will be increased, AS on the fig. Id will lose stability and 
become the same skyrmion as the others around. 




Figure 2: Calculated structure of pure AS. (a) - in zero 
field; (b) - in the field h = 0.4. 

For studying the structure of undistorted AS core we 
used the formula (7) for fixing boundary conditions. Cal- 



3 



ciliated AS structure is presented on the fig. 2. It is notice- 
able that when the external field increases, then period 
increases too and areas, which are magnetized against 
the field become thinner. Dependence on the fig. 3 shows 
that (10) is suitable even for the case of only two turns. 
Values for k = 3.5 does not match a little to the de- 
pendency, but it was made for being closer to the limit 
of accuracy of the computational approach. Required 
accuracy grows in proportion to k 2 and for k > 3 is al- 
ready comparable with calculating accuracy. By using 
a straight line approximation (fig. 3), according to (10), 
determining appropriate value for k = 1, it is possible to 
calculate the core energy. The obtained dependence 
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in common with the formula (10) may be used for high 
accurate calculation of AS energy with any quantity of 
turns. 
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Figure 3: AS energy dependence (relative to the helicoid) 
from the number of turns, in logarithmic scale. Thin line 
- in zero field, thick line - in the field h = 0.4. 



New equilibrium state - hexagonal lattice of 
Archimedean spirals (ASL) with different number 
of turns was found using the method of average en- 
ergy density minimization with periodical boundary 
conditions. Such structures remain stable in small 
perturbations and do not transform into helicoid. Super- 
position of identical precomputed solitary AS was chosen 
for the initial state in such a way that four centres were 
set in the corners of rectangular cell (fig. 4) and one - in 
the middle. In the process of calculation turns changed 
shape to hexagonal and the cores of spirals distorted and 
slightly elongated. The central core rotated noticeably 



relative to the corner ones. On the fig. 4 a structure of 
equilibrium lattice with five turns is shown. Though this 
phase is similar to SD lattices, which were observed in a 
certain kind of ferrimagnetic films [16, 26], its existence 
is due to the local Dzyaloshinskii-Moriya interaction but 
not to the long-ranging magnetic dipole forces. 

We made a numerical analysis of the energy of different 
two-dimensional structure with zero temperature and it 
showed the following. In zero magnetic field the average 
density of energy (in units Eq/L 2 d ) of Dzyaloshinskii he- 
licoidal phase gb = —19.74, for ASL - gas — —19.58, 
for skyrmion lattice - as = —18.31. Negative sign 
in expression is tied with negative sign of energy of 
Dzyaloshinskii-Moriya interaction. In the magnetic field 
when $J h < 0.2 values of these energies are ordered as 
o~d < & as < &Si ASL loses its stability when h « 0.2. 
That is why ASL may be a basic state only when tem- 
peratures are different from zero, in the range of fields 
^ ft < 0.2. As the number of ASL degrees of free- 
dom is much higher than the number of one-dimensional 
helicoid degrees of freedom then its free energy may be 
smaller. For the verification of this statement we should 
calculate system free energy (with intermediary temper- 
atures where AS are actually observed [3]). At present, 
these calculations cannot be fully completed neither the- 
oretically nor numerically. ASL as a metastable state 
may be obtained in periodical magnetic field as it is used 
for creation of SD lattices [26]. 




Figure 4: Calculated ASL structure with five turns in 
zero field. Selected domain is the computational grid 
area, which is the unit cell. 
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Figure 5: Color scheme, which determines the direction 
of the magnetization vector in three dimensions. 
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